Dynamical decoupling and dephasing in interacting two-level systems 
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We implement dynamical decoupling techniques to mitigate noise and enhance the lifetime of an 
entangled state that is formed in a superconducting flux qubit coupled to a microscopic two-level 
system. By rapidly changing the qubit 's transition frequency relative to the two-level system, we 
realize a refocusing pulse that reduces dephasing due to fluctuations in the transition frequencies, 
thereby improving the coherence time of the entangled state. The coupling coherence is further 
enhanced when applying multiple refocusing pulses, in agreement with our 1// noise model. The 
results are applicable to any two-qubit system with transverse coupling, and they highlight the 
potential of decoupling techniques for improving two-qubit gate fidelities, an essential prerequisite 
for implementing fault-tolerant quantum computing. 



A universal set of quantum gates, sufficient for imple- 
menting any quantum algorithm, consists of a two-qubit 
entangling gate together with single-qubit rotations [1 . 
However, fault-tolerant quantum computing with error- 
correcting protocols sets strict limits on the allowable 
error rate of each gate. Initial work focused on per- 
fecting single-qubit gates [2H1], but in recent years there 
has been progress on characterizing two-qubit gate op- 
erations [5^. In superconducting systems, two-qubit 
gates have been implemented in a variety of ways, for 
example through geometric couplings [8, 9 , tunable cou- 
pling elements [lOlITT], microwave resonators tUlIIS], or 
with microwave- induced interactions [6\ [MUTT] . Regard- 
less of the nature of the coupling, any variations of the 
qubit frequencies or in the coupling parameter during the 
two-qubit interaction leads to dephasing of the entangled 
state, and puts an upper limit on the obtainable gate fi- 
delity [5 . 

For single qubit s, dephasing due to low- frequency fluc- 
tuations in the precession frequency is routinely reduced 
with refocusing techniques [18-20^. originally developed 
in nuclear magnetic resonance [21]. In this work, we ap- 
ply similar techniques to improve the coherence of an 
entangled state formed between a flux qubit and a mi- 
croscopic two-level system (TLS). The refocusing pulse 
is implemented by rapidly changing the qubit frequency 
relative to the TLS, thereby acquiring a phase shift [22 . 
When the phase shift equals tt, the pulse refocuses the 
incoherent evolution of the coupled qubit-TLS system, 
giving a fourfold improvement of its coherence time. 
We further prolong the decay times by applying multi- 
ple refocusing pulses [23, 24^ , thus extending dynamical- 
decoupling techniques [25] to correct for dephasing of en- 
tangled states. The results are first steps towards imple- 
menting error-correcting composite gate pulses [26, 27] 
and optimal control methods [28 , schemes with strong 
potential for improving two-qubit gate operations. 
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FIG. 1: (a) Spectroscopy of the qubit-TLS system. The qubit 
and TLS are resonant at / = 7.08 GHz, where the spectrum 
has an anticrossing with splitting S — 76 MHz. The inset 
shows the qubit spectrum over a larger range, with the red 
circle indicating the region of interest, (b) Pulse sequence for 
probing the qubit-TLS interactions. The vr-pulse generates 
a qubit excitation, which is coherently exchanged back and 
forth between qubit and TLS during the interaction time n. 
(c) Coherent oscillations between qubit and TLS, measured 
using the pulse sequence shown in (b). High switching prob- 
ability Psw corresponds to the qubit's ground state |0), low 
Psw to the qubit's excited state (d-e) Characteristic de- 
cay time tdecay and oscillation frequency /osc, extracted from 
the data in (c). The oscillations decay faster for 6^ ^ 0, 
a consequence of the increased sensitivity dfosc/d^ to flux 
noise. 



We use a flux qubit [29] , consisting of a superconduct- 
ing loop interrupted by four Josephson junctions (see 
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Ref. [24] for a detailed description of the device). The 
qubit 's diabatic states correspond to clockwise and coun- 
terclockwise persistent currents ±/p, with Ip = 180 nA. 
The inset of Fig.jlja) shows a spectrum of the device ver- 
sus external flux, with <l>qb defined as = ^-\-^o/2 and 
$0 = h/2e. The qubit frequency follows /qb = VA^ + e^, 
where the tunnel coupling A = 5.4 GHz is set by the de- 
sign parameters and the energy detuning £ = 2/p<l>qb//i 
is controlled by the applied fiux. The device is embedded 
in a SQUID, which is used as a sensitive magnetometer 
for qubit read-out [18 . 

At <l>qb = <!>* = ±4.15 m<l>o, the qubit becomes resonant 
with a TLS |30^. The microscopic nature of the TLS 
is unknown, but studies of two-level systems in sim- 
ilar qubit designs show that the most likely origin is 
an electric dipole in one of the tunnel junctions [31 . 
Figure [TJa) shows a magnification of the region around 
— 4.15m^0 7 revealing a clear anticrossing with splitting 
S = 76 MHz. We describe the system using the four 
states {|0^), |1^), |0e), |le)}, where (0, 1) are the qubit 
energy eigenstates and (^, e) refer to the ground and ex- 
cited state of the TLS. On resonance, \lg) and |0e) are 
degenerate and coupled by the coupling energy hS. To 
characterize the coupling, we use the pulse scheme de- 
picted in Fig.jljb) [32 . Starting with both qubit and 
TLS in their ground states (|0^)), we rapidly shift the 
fiux to a position = <l>qb = 1.2m^o where the 

qubit frequency /qb = 6.3 GHz is far detuned from the 
TLS. By applying a microwave pulse, resonant with /qb, 
we perform a 7r-rotation on the qubit and put the system 
in \lg). We then rapidly shift to a value close to zero, 
effectively turning on the interaction whereupon the 
system will oscillate between \lg) and |0e). After a time 
Ti , the interaction is turned off by shifting 5^ away from 
zero, and we measure the final qubit state by applying 
a read-out pulse to the SQUID. Since the measurement 
outcome is stochastic, we repeat the sequence a few thou- 
sand times to acquire sufficient statistics to estimate the 
SQUID switching probability Psw and thereby the qubit 
state. 

Figure [TJc) shows the qubit state after the pulse se- 
quence, measured versus interaction time ri and ffux de- 
tuning S^. At = and for n = 1/(25') = 6 ns, 
the pulse sequence implements an iSWAP gate between 
qubit and TLS, taking |0e) ^ i|l^) and |1^) ^ i\Oe) [32]- 
[34]. The characteristic decay time of the oscillations is 
shown in Fig.jljd). The oscillations persist the longest 
at = 0; at this point, the decay time is ~ 800 ns. 
However, as is moved away from the optimal point, 
the decay time quickly decreases to zero. We attribute 
the reduction in coupling coherence to low-frequency ffux 
noise, present in all superconducting devices [35]. When 
7^ 0, ffuctuations in induce variations in the ef- 
fective coupling frequency /osc pig.jlje)], leading to de- 
phasing of the entangled state. 

For single qubit s, dephasing due to low- frequency ffuc- 
tuations of the qubit frequency can be reduced in a Hahn- 
echo experiment [21 . By applying a 7r-pulse after a time 



t of dephasing, the qubit 's noise-induced evolution will 
reverse directions and refocus at time 2t, provided that 
the ffuctuations are slow on the time scale 2t |24l [36] . 
Here, our goal is to extend such single-qubit refocusing 
techniques to the mitigation of noise in coupled systems 
with multiple qubits, which requires implementing refo- 
cusing pulses for entangled states. Note that the purpose 
here is to increase coherence times, as opposed to turning 
off unwanted couplings [371 131] • 

We start by describing the system's dynamics. Fol- 
lowing Refs. [30l[3TJ|35, we write the total Hamiltonian 

as ^ = i^qb + ^TLS + ^int, with i^qb = -(/l/2)/qb , 

^TLS = — (^/2)/tls <^J^^, and with the interactions de- 
scribed by ^int = -{h/2)Safa^^^. Here, a^J^^ are Pauh 
operators for the qubit, aj^^ are TLS operators and /tls 
is the TLS frequency. To focus on the interactions be- 
tween the qubit and the TLS, we restrict the discussion 
to the subspace spanned by the states {|1^), |0e)}. In the 
rotating frame of the TLS, the subspace Hamiltonian be- 
comes 

Hsub = -'^{SfaT'' + SaT'^), (1) 

where 6f = /tls — /qb and are subspace Pauli op- 
erators. The dynamics of Eq. (fTl) can be visualized on 
a Bloch sphere, with the north and south poles corre- 
sponding to \lg) and |0e), respectively, and with S and 
df representing the length of torque vectors along the x- 
and 2:- axes [see Fig.[2|b)]. The frequency of the coherent 
oscillations seen in Fig.mc) is then given by the effective 
coupling strength 

/osc = VW^. (2) 

which is plotted together with the data in Fig.jlje). 

With the coupling dynamics described by Eq. ([T]), we 
discuss the details of the refocusing sequence, shown in 
Fig.[2|a-b). The system is brought into the {|1^), |0e)} 
subspace by applying a 7r-pulse to the qubit [step I in 
Figs.[2|a-b)], followed by a non-adiabatic shift in to 
bring the qubit and TLS close to resonance. \lg) is not 
an eigenstate of the coupled system, so the interaction S 
will cause the system to rotate around the x-axis, oscil- 
lating between \lg) and |0e). Low- frequency ffuctuations 
in the effective coupling strength will cause the Bloch 
state vector to fan out (over many realizations of the ex- 
periment), and the system loses its phase coherence (step 

II). 

The refocusing pulse is now implemented by apply- 
ing a ffux shift pulse that rapidly detunes the qubit and 
the TLS to Sf = 550MHz. With \Sf \ > \S\, the state 
vector is effectively rotating around the z-axis (step HI) 
[19], and we realize a tt rotation by setting the pulse 
duration Trefocus = 0.5/^/. The system is then rapidly 
brought back into resonance (step IV), and the state 
vector continues to rotate around the x-axis. The in- 
homogeneous broadening that caused the state vector to 
diffuse during the ffrst interval ri will now realign them 
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FIG. 2: (a) Pulse sequence and (b) Bloch sphere repre- 
sentation of the refocusing protocol. The blue arrows are 
state vectors, while the red arrows represent the Hamilto- 
nian in Eq. (fTl). After the qubit tt pulse, the system enters 
the {|1^), \0e)} subspace (step I). The coupling S rotates the 
state vector around the x-axis, but due to noise in the effective 
coupling the state vector fans out (II). The rapid flux pulse 
^refocus gcncratcs a large frequency detuning Sf, the system 
will start rotating around the z-axis (III) and eventually com- 
plete a TT-rotation (IV). The inhomogeneous broadening now 
refocuses the state vector, giving an echo at V. (c-d) Evolu- 
tion of the qubit-TLS system, measured with and without a 
refocusing pulse. A clear echo appears after the refocusing 
pulse, with a maximum close to T2 = n. The traces were 
taken at 6^ = -72 fi^o. 



again. The refocusing is complete after a time r2 = ri 
(step V) . Figures |2|c-d) illustrate the result of the refo- 
cusing sequence. Without the refocusing pulse [Fig.[2|c)], 
the coherent oscillations between \lg) and |0e) decay al- 
most completely after 100 ns. When inserting a refocus- 
ing pulse at Ti = 100 ns, the oscillations start to revive, 
eventually forming an echo at r2 = ri. 

The revival of phase coherence seen in Fig.[2|d) re- 
quires careful calibration of the refocusing pulse. FigurejS] 
shows an example of a calibration experiment, where 
we fix Ti = 97 ns and 5f = 550 MHz and measure re- 
focused oscillations versus the refocusing time Trefocus- 
The data shows strong oscillations whenever the refocus- 
ing pulse rotates the state vector by an odd integer of tt, 
i.e. when Trefocus = (2n + 1) x 0.5/6f^ in agreement with 
the schematics discussed in Fig.[2[b). 

We now turn to investigating the decoherence mecha- 
nisms and determining the performance of the refocusing 
protocol. For Gaussian-distributed dephasing noise, we 
expect the amplitude h{t) of the coherent oscillations to 
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FIG. 3: Calibration of the refocusing pulse, measured by 
fixing n = 97 ns and looking for the echo signal by sweep- 
ing T2. The echo appears every time the refocusing pulse 
generates a rotation by an odd integer of tt. The first five 
refocusing conditions occur at Trefocus = (2n + 1) x 0.5/ 5f = 
0.9, 2.7, 4.5, 6.4, and 8.2 ns, with Sf = 550 MHz. 



decay as [39] 



h{t) = exp[-t/fi] exp[-(t/T^,N)']. 



(3) 



The exponential decay constant Ti is due to energy relax- 
ation, while T(^^N represents the dephasing with N refo- 
cusing pulses. At = 0, we measure a pure exponential 
decay with time constant Ti = 800 ns, which is shorter 
than the relaxation time of both the qubit {T^^ = 10 /is) 
and the TLS {T^^^ = 1/is). However, to get an ex- 
pression for Ti, we need to consider all the possible ab- 
sorption/emission rates in the full four- level system [40] , 
In the relevant situation hS <C /cbT <C hfTLS, ^/qb, we 
have 



l/Ti 



(r± + r^), (4) 



where F-^ (r^) represents relaxation (excitation) be- 
tween the two energy eigenstates |±) = (|0^) ± \ le))/V2 
of Eq. ([1]), with energy splitting /i/osc- The polarization 
rate T± + F^ = 5'^(/osc)/2 depends on the noise power 
S± that couples transversely to the diagonalized subspace 
Hamiltonian in Eq. (IT]), which for 6^ = corresponds to 
fluctuations ^'^/(/osc) frequency detuning Sf [20\ . 

Using Eq. ^ and the measured values of Ti , T^^ and 
Ti^LS, we get %(/ = 76 MHz) = 2.8 x lOSad/s. We 
can not distinguish whether this noise comes from fluc- 
tuations in /qb, /tls or a combination thereof, but we 
note that the measured value is a few times larger than 
fluctuations in /qb expected from flux noise. From inde- 
pendent measurements of the flux noise power , in 
the same device, we have Sf^^ = 5'<^>qb(^/qb/S^qb)^ = 
1.1 X lO^rad/s at / = 76 MHz and <l>qb = -4.15m<l>o 

Away from 6^ = 0, the decay envelope becomes Gaus- 
sian, and we extract the dephasing time T^^n by fitting 
the data to Eq. ([3|, assuming a constant relaxation time 
Ti = 800 ns. The extracted decay times versus flux 6^ 
are shown in Fig.|4]^a). The refocusing sequence gives 
considerably longer decay times over the full range of the 
measurement except around ^<l> = 0, where the decay is 



4 




(5$ = -48 /^$c 
-•- 5$ = -60 /x$c 
-•- 5$ = -72 

= -83 



-100 -50 

Flux detuning [m^oI 



12345678 
# of refocusing pulses N 



T 1 1 1 1 1 r 




No refocusing (A/=0) 



>^ Refocusing pulses . 



100 200 300 400 500 600 700 
Time [ns] 

FIG. 4: (a) Decay times of the coherent oscillations, mea- 
sured with (A^=l) and without (A^=0) a refocusing pulse. At 
6^ = 0, the decay is limited by energy relaxation, but the 
coherence times decrease away from = due to increased 
sensitivity to flux noise. For large |^$|, the refocusing se- 
quence improves the decay time by more than a factor of four. 
The inset show examples of decay envelope /i(t), measured at 
— -60/i$o- For the A^=l case, we use n = T2 = t/2. 
(b) Decay times of multi-pulse refocusing sequences, showing 
the improvement coherence as the number of pulses increases, 
(c-d) Time evolution of the qubit-TLS system, measured with 
and without refocusing pulses. Note the echo signals appear- 
ing after each refocusing pulse, giving a strong enhancement 
of the coherence time. The data is taken at 6^ — —84 /i$o- 



limited by Ti. Note that to capture both the exponen- 
tial and the Gaussian decay, we plot the time Tg for the 
envelope to decrease by a factor 1/e. Examples of de- 
cay envelopes together with fits are shown in the inset of 
Fig.|4ja), measured with and without a refocusing pulse 
at 6^ = 

We model the decreased phase coherence away from 
J<l> = in terms of flux noise. In analogy with co- 
herence measurements on single flux qubits we as- 
sume l//-type fluctuations in ^, with noise spectrum 
S^{uj) = A^/\uj\. The fiux noise couples to the oscillation 
frequency /osc through Eq. ([2|, leading to the dephasing 
rate 



1/T^,N = 27ry^CNA^\dfosc/d^\. 



(5) 



Here, cn=o = ln(l/cJiow^) and cn=i = ln(2) relate to the 
filtering properties of the pulse sequence [20l[2j, with the 
low-frequency cut-off a;iow/27r = 1 Hz fixed by the mea- 
surement protocol. The solid lines in Fig.[4]^a) are fits to 



Eq. ([5|, using a single fitting parameter A^ = (1.4 /i^o)^- 
This amount of flux noise is consistent with previous re- 
sults miig. 

The overall good agreement between Eq. ([5| and the 
data verifies the noise model and further confirms the va- 
lidity of the refocusing sequence. However, for the range 
5^ > — 30/i<l>o, the refocused data shows slightly lower 
coherence times than expected from the model. We at- 
tribute this to the finite rise time of our shift pulses. The 
refocusing sequence requires the frequency sweep rate 
df/dt to be fast compared to the interaction timescale 
l/^^lOns (to make the shifts non-adiabatic), but slow 
compared to the qubit precession time l//qb ^ 0.2 ns 
(to avoid driving the system out of the {|1^), |0e)} sub- 
space). The constraints can be phrased in terms of the 
probability of undergoing Landau- Zener transitions, giv- 
ing < df/dt < Using pulses with 1.5ns 
Gaussian rise time, we have df/dt ~ 550 MHz/1.5 ns = 
(606 MHz)^, and on average the constraints are well ful- 
filled. However, df/dt is lower during the slowest parts 
of the shift pulse (the beginning and the end), and ar- 
tifacts due to imperfect non-adiabaticity appear when 
these parts of the pulse occur where the effective cou- 
pling is the strongest (at 5^ = 0). The limited non- 
adiabaticity is also the reason for the slight asymmetry 
around (5<l> = in Fig.jljc) [22] . 

We now extend the refocusing technique to imple- 
ment dynamical decoupling protocols with multi-pulse 
sequences. For 1 / /-type noise, it has been shown that the 
Carr-Purcell sequence |23], consisting of equally spaced 
TT-rotations, improves coherence times by filtering the 
noise at low frequencies [2T, ^361 HH [45]. Figures |4jc- 
d) show the coherent evolution of the system when re- 
peatedly applying refocusing pulses. Echo signals form 
between each pair of tt pulses, giving considerable longer 
coherence times compared to the N = case. The in- 
crease in decay time with the number of refocusing pulses 
N is plotted in Fig.Qb), measured for a few different 
values of 6^. The improvement is consistent with the fil- 
tering properties of the pulse sequence; for larger A/", the 
filter cut-off frequency increases and, since the noise is of 
l//-type, the total noise power leading to dephasing is 
reduced. 

To summarize, we have implemented refocusing and 
dynamical decoupling techniques to correct for noise and 
improve the lifetime of entangled two-level systems. Al- 
though implemented between a flux qubit and a mi- 
croscopic two-level system, the method applies to any 
transversely coupled spin- 1/2 systems where the relative 
frequency detuning can be controlled. We expect the 
findings to be of importance when developing decoupling 
techniques to improve two-qubit gate fidelities. 

We thank K. Harrabi for assistance with device fabri- 
cation and O. Zwier, X. Jin and E. Paladino for helpful 
discussions. This work was supported by NICT. 
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